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The controlled switching between two quasi-stable Neel states in adsorbed anti-ferromagnetic Fe 
chains has recently been achieved by Loth et al [Science 335, 196 (2012)]. In order to rationalize 
their data, we evaluate the rate of tunneling electron-induced switching between the Neel states. 
Good agreement is found with the experiment permitting us to identify three switching mechanisms: 
(i) low-bias direct electron-induced transitions, (ii) intermediate-bias switching via spin-wave-like 
excitation, and (iii) high-bias transitions mediated by domain wall formation. Spin correlations 
in the anti-ferromagnetic chains are the switching driving force leading to a marked chain-size 
dependence. 

PACS numbers: 68.37.Ef, 72.10.-d, 73.23.-b, 72.25.-b 



The search for nano-scale electronic devices has 
prompted intense research in the field of nano- 
magnetism. The use of spin as the information convey- 
ing entity has steered much excitement due to its ex- 
traordinary properties of information storage, speed and 
low-energy consumption [3-[3] ■ Miniaturization is quickly 
proceeding, reaching very small domain- wall devices [4|, 
atomic-size devices [H| and the realm of molecular de- 
vices. [1,04^ Among all these possibilities, antiferromag- 
netically (AFM) coupled devices have recently received 
a lot of attention. The AFM characteristics make these 
devices very well fitted for quantum com put ation since 
they naturally involved entangled states. [lCt [llT ] More- 
over, the storage in AFM devices is particularly robust 
due to the lack of a total magnetic moment. However, 
this robustness has deterred their use because changing 
their magnetic state becomes difficult 12 1. 



Recently, Loth and co-workers succeeded in control- 
lably switching the spin states of AFM atomic chains. [12 1 
Two quasi-stable Neel states, exhibiting alternating spin 
directions on the atoms along the chain, were evidenced 
in Fe chains adsorbed on a CuN/Cu(100) surface. Loth 
and co-workers showed that the Neel states can be 
switched by tunneling electrons injected from a polarized 
Scanning Tunneling Microscope (STM) tip into one of the 
atoms of the chain. This demonstrated the possibility 
of storing information on atomic scale anti-ferromagnet. 
Theoretical predictions show that writing and reading 
spin states entail fundamental problems associated with 
the quantum nature of the process. [HI Spin manipu- 
lation by tunneling electrons has been pictured as due 



to a spin-torque mechanism where spin angular momen- 
tum from the electron is transferred into the atomic spin 
system. [I4l4l6| However, due to the lack of magnetic mo- 
ment in AFM systems, spin manipulation must follow a 
different mechanism. Unavoidably, spin manipulations 



and excitations are closely related. [17| Indeed, switch- 
ing between Neel states has been experimentally associ- 
ated with overcoming an activation energy. [l2| However, 
in the case of degenerate Neel states, resonant transi- 
tions should be expected. Hence, the experimental data 
raise many questions regarding the possibility of resonant 
switching, the efficiency of activated switching, the na- 
ture of the involved excitations and the physics at play 
in AFM spin torque. In summary, a complete view of the 
switching process is missing. 

In this Letter, we reveal the switching mechanisms at 
play in the experiment of Rcf. 



121. The mechanisms 



turn out to be rich and closely related with the excita- 
tion spectra of the AFM chain. Their understanding give 
us a handle on the parameters controlling the switching 
process. Here, we show that the correlated spin nature 
of AFM quantal chains is at the origin of the transition 
between Neel states, and from this, we deduce the behav- 
ior of the switching rate with respect to applied bias and 
chain size. Hence, our theory shows that the ability of 
switching states are intrinsic to AFM-correlated atomic 
systems. 

One Fe atom on CuN/Cu(100) is characterized by an 
5 = 2 spin 18, 19| with a large magnetic anisotropv jlj| . 
The easy axis lies along the line of N-atoms of the surface. 
Experimental 12 , 2(| and theoretical 21 1 evidence show 
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that in a first approximation, a chain of transition metal 
atoms on CuN/Cu(100) is an ensemble of weakly inter- 
acting atoms such that the Fe chain can be described as a 
set of S = 2 spins with magnetic anisotropy and coupled 
by Heisenberg exchange. This is partially due to the de- 
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coupling properties of the CuN layer [14|,[16|,[20| and to 
the considerable distance between Fe atoms. Hence, the 
system can be described by the magnetic Hamiltonian: 
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+ ]T (g^ B B ■ Si + DSl + E{Sl x - £!?„)) , (1) 
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where Si is the spin of the atom i (i = 1, N) and Si lU its 
projection on the u-axis. D and E are the longitudinal 
and transversal anisotropy coefficients (D < 0). B is a 
macroscopic magnetic field applied to the system, along 
the z-axis, similarly to the experiment [T3 | . 

Neel states are broken-symmetry solutions of Hamilto- 
nian (Q}. Hence, Hamiltonian (p} cannot represent the 
Fe chains of Ref. 12j in the absence of an inhomogcnc- 
ity that breaks the symmetry. In the present study, we 
introduced a phenomenological term to enforce the Neel 
magnetic structure. Hence, the magnetic Hamiltonian 
becomes: 
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FIG. 1: (a) Setup used in the experiment of Ref. [12[ ]. 
A few atoms of an Fe (large gray balls) monoatomic chain 
on CuN/Cu(100) are shown under an STM tip that injects 
spin-polarized electrons. The Fe-N bonds are also depicted, 
the presented arrangement corresponds to the atomic coordi- 
nates of a DFT calculation. The chain spins alternate along 
the chain axis to form a Neel state, (b) Efficient excitation 
and de-excitation schemes leading to the switching between 
Neel-states. D refers to the direct process between the quasi- 
degenerate Neel-states and / refers to the indirect one that 
involves excitation and de-excitation mechanisms. The up ar- 
row on the left shows the excitation by a tunneling electron 
whereas the down arrows on the right show the relaxation in- 
duced by substrate electrons. All possible de-excitation cas- 
cades are included, signaled here by a few de-excitation paths. 
The depicted spectrum corresponds to a chain of 10 Fe atoms. 
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A small inhomogeneous field, Bi n h, of 0.1 T acting 
on one of the atoms of the chain (here an end atom) 
is enough to split the ground states into two Necl-likc 
states (noted 1 and 2) that contain contributions from 
many spin configurations. This added perturbation is in- 
deed small since the two Neel states are only ~ 50/xeV 
away from each other. The inhomogeneous term can be 
thought to be representative of various effects: small in- 
homogeneities of the surface, an inhomogeneity of the ex- 
ternal B field, dephasing effects or a field induced by the 
polarized tip of the STM [Hj]. Following Loth ct al. [H, 
the parameters in Hamiltonian (Q]) were partially deter- 
mined by fitting the experimental Fc2 excitation energies. 
Hence, similarly to the fitting of the Ising-model param- 
eters in Ref. [12j, we set E = and obtained J = 1.6 
meV and D = —1.34 meV (23[. Here we consider even- 
numbered chains and the Hamiltonian ([2]) was diagonal- 
ized in a basis formed by direct products of local spin 
states (spin configurations): nfc|M(fc)), where \M(k)) is 
an eigenstate of the projection on the z-axis of the lo- 
cal spin at atom k, Sk,z- For the longer chains, diago- 
nalization methods specific to sparse matrices had to be 
used 0. 

We describe the magnetic excitation during tunneling 
using the strong coupling approach of Ref. [19( in which 
the electron strongly interacts with the atom it tunnels 



through. Let us consider the system schematized in Fig.Q] 
(a): the STM tip is standing above the first atom of a 
chain in one of the two Neel states and the injected elec- 
tron is tunneling through it. Within the sudden approx- 
imation the electron transmission amplitude operator is 
equal to: [27 1 
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where St is the spin of the compound system formed by 
the tunneling electron and the corresponding Fe atom. 
^Tip^rSub ^ S ^ ne electron transmission amplitude in the 
St symmetry from tip to substrate. Here we use St — 9 
revealed as the dominant tunneling channel by density 
functional theory studies in Ref. [19| for an isolated Fe 
atom on CuN/Cu(100). The scattering amplitude be- 
tween the chain magnetic states is then obtained as 
the matrix element of the amplitude ([3]) between initial, 
\i) = \oi)\4>i) and final states, |/} = \<Jf)\4>f), written as 
direct products of the scattering electron state of spin a 
and the chain state \4>), eigenstate of Hamiltonian [2] The 
probability, P(i — > /), for a transition from the initial 
state |i) to the final state |/) associated to the tunnel- 
ing of an electron from tip to substrate is proportional 
to \T^ S J 
5/2,M T |i)| 2 . 



and to \E MT (f\ S T = 5/2,M T )(S T = 
2 factor corresponds to a 
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FIG. 2: Tunneling electron-induced switching rate between 
Neel states in s _1 as a function of the applied bias in mV 
for Fe chains of different length. The higher rates correspond 
to the Fe dimer, and the rate decreases as the chain length 
increases. For comparison the experimental results p2] are 
plotted for transitions from low-current to high (dots) and 
high to low (diamonds). 



global tunneling probability, whereas the second factor 
yields the relative importance of the inelastic channels 
in the tunneling process. The (i— >f) probability is very 
large if there is a strong overlap between the intermedi- 
ate state of St tunneling symmetry and both the initial 
and final states. 

The electron can only induce a change of one unit in 
the spin projection of the atom it tunnels through. As 
a consequence, in a zero-order view of the Neel states, 
described as chains of atoms with alternating spins, 
electron-induced switching between pure Neel states does 
not exist. However, AFM chains described with Hciscn- 
berg couplings are strongly correlated 25, [261, i.e. the 



two Neel- like states, eigenstates of Hamiltonian ([2]), con- 
tain small components over a very large number of differ- 
ent configurations of atomic spins. Direct, quasi-resonant 
electron-induced transitions between Neel-like states are 
then possible. These quantal transitions that do not in- 
volve spin-flip of the tunneling electron are mediated by 
correlation. However, since they involve the small com- 
ponents in the state expansion over spin configurations, 
their probability is weak. Furthermore, if the length of 
the chain is increased, the direct 1— s-2 transition probabil- 
ity decreases rapidly. This is clearly seen in Fig. [2] which 
shows the rate for the 1— >2 transition induced by tunnel- 
ing electrons as a function of the STM bias. At low bias 
(below 6 mV), only direct 1— >-2 transitions are possible 
and lead to a very small transition probability for long 
chains that very rapidly decreases as the chain length in- 
creases. Typically, the direct 1— >2 excitation probability 
per electron is equal to 0.9 x 10" 7 (1.5 x 10~ n ) for Fe6 
(Fe 8 ). 



Above 6 mV, the transition rate increases drastically 
due to indirect transitions: the system initially in state 1 
is excited by collision with the tunneling electron into 
a higher- lying state i, and state i later decays into 
state 2 with a finite probability (Fig. [1] (6)). The de- 
cay process takes place via electron-hole pair excita- 
tion of the substrate and it is the inverse of the exci- 
tation process discussed above. [l6| The above strong- 
coupling approach successfully accounts for excited spin 
state lifetimes. Hence, we can complete the full 

cxcitation/dc-excitation dynamics, by considering now 

ec. 



the de-excitation probability P® ec 
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where P{i — > f) has been defined earlier and AEij = 
Ei — Ef is the decay energy. Since there is no long-lived 
excited state in the system (except state 2), an excita- 
tion of state i, followed by successive decays according to 
dU, eventually results in the population of the two Neel 
states, iV,(l) and Ni(2). The 1— >2 transition rate for a 
bias V is then given by: 

Rate{\^2) = C[(eV - AE 2 i)P(l ->• 2) 

+ J2( eV ~ ->• 0W)]. (5) 

i>2 

The C coefficient is the normalization constant for the 
total flux of electrons in a given experiment [29[ which is 
bias independent for small bias. The first term in Eq. (|5|) 
corresponds to the direct transitions already discussed 
and the sum over i to the indirect transitions. 

Figure [2] shows the computed rates for Fe„ chains 
(n = 2 — 10) compared with the experimental results for 
Fes- [IH The theoretical transition rate decreases with 
chain length and from there one can conclude that di- 
rect transitions are practically impossible for very long 
chains whereas the indirect process should be accessi- 
ble in a broad range of lengths. Our results on Fig. [5] 
compare very well with the experimental data for Feg, 
also shown as dots and diamonds. Dots in Fig. [2] show 
the experimental transitions rates from low current to 
high current, and diamonds from high to low currents. 
This corresponds to 1 —5- 2 transitions and 2 — >• 1 respec- 
tively. The experimental results in Fig. [5] were obtained 
with a polarized tip, whereas the theoretical results were 
obtained with a non-polarized tip. Polarizing the tip 
influences the transition rates, due to the intermediate 
St=5/2 symmetry involved in the tunneling. Inelastic 
transitions are favored when the tip and the atom under 
the tip have opposite spin directions; the ratio between 
the theoretical switching rates is typically around 3 for a 
tip polarization equal to 0.3 which explains the asymme- 
try between dots and diamonds in Fig. [2 

The theoretical switching rates show an abrupt change 
at <~ 12 meV. This change is more clearly seen in Fes 
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and Feio chains. By studying the spectra of excitations 
for these chains, we can clearly separate two distinct re- 
gions in the rate due to two sets of excitations different 
in nature. The first region, between 6 and 12 mcV, cor- 
responds to transitions where M changes in one of the 
end atoms of the chain; these are very efficient in shorter 
chains such as Fe6 . This type of excitation is a quantized 
spinwave of the finite chain (26^ . The second region, be- 
yond 12 meV, corresponds to transitions that mix several 
configurations with two opposite anti-ferromagnetic do- 
mains, they lead then to domain-wall formation. As the 
chain length increases, domain- wall mediated transitions 
become more important than spinwave-mediatcd transi- 
tions, explaining the clear upturn beyond 12 meV for Fes 
and Feio- There is also a dependence of the transition 
rate on the position of the tip along the chain. Indeed, 
this position effect mixes with the polarization effect due 
to the alternating spin directions along the chain in the 
Neel state as well as with the effect of the position of 
the inhomogeneous term (2) along the chain. Besides 
the interplay with the polarization effects, end atoms are 
more efficient for moderate bias and central atoms for 
high voltages. 

The transition path 1 — > i — > 2 for Fes is analyzed 
in Fig. [3] The different contributions to the switching 
rate, Eq. ((5|), are presented as function of the excita- 
tion energy: the primary excitation probability P(l — > i) 
(black circles), the branching ratio toward state 2 ( A^(2), 
red line) and the total indirect excitation probability 
P(l — > i)Ni(2), (green diamonds). Many states i are ex- 
cited by a tunneling electron, though many of them with 
a small probability that rou ghly exponentially decreases 
with the excitation energy [26J. The largest excitation 
probability corresponds to states i where the spin of the 
atom under STM tip changes by AM rj ±1. All the 
other states are excited via correlation, i.e. the fact that 
states 1 and i are not associated to a single configuration 
of local spins [26| but to a mixing of a large number of 
them from a configuration interactions point of view. As 
for the decay of the excited states, i, the lower-energy 
states decay preferentially toward Neel state 1 or 2 de- 
pending on which state they are configurationally closer. 
For high-lying states, correlation becomes stronger and 
the excited states decay roughly equally to the two Neel 
states. Domain walls in particular decay roughly equally 
to both Neel states. As a conclusion from Fig. [3l the 
global transition 1 — > i — > 2 requires that both excitation 
and de-excitation are sizeable, i.e. requires balancing the 
1 — > i and i — > 2 probabilities. 

The relative weight of the various switching processes 
vary with chain length and this can be rationalized con- 
sidering that the distance between the two Neel states 
in terms of spin-configuration changes is increasing with 
chain length. For the direct quasi-resonant process at 
small bias, one needs a strong configuration mixing be- 
tween the two Neel states and this quickly decreases with 



10" 



10" 



o 
- 



10" 



10" 




10 15 

Energy (meV) 



20 



FIG. 3: Probability as a function of the i-level energy for 
the excitation from the Neel state 1 to level i: P(l — > i) 
(black circles); branching ratio toward state 2 once level i is 
populated: Ni(2) (red line); and total indirect switching prob- 
ability between Neel state 1 and 2 via an intermediate i-level: 
P(l —¥ i)Ni(2) (green diamonds). Many states i are excited 
by an electron, with a probability roughly exponentially de- 
creasing with the excitation energy, but the final switching 
probability is strongly modulated by the de-excitation prob- 
ability (branching ratio). 



chain length. The process around 6-8 meV is associated 
to AM « ±1 transitions in the atom just under the tip; 
the corresponding excited state is still very close to the 
initial Neel state and therefore the decay to the other Neel 
state is difficult and rapidly decreases with chain length 
increase. The process around 12-13 meV is associated 
to domain wall formation and behaves differently. First, 
several states contribute and their number increases with 
chain length; second, their decay equally populates the 
two Neel states. So, even, if their excitation probability 
from the initial Neel state is not very high, they succeed 
in dominating the indirect process for long chains. 

In summary, our calculations show three different 
regimes in the tunnelin g el ectron induced switching of 
the AFM chains in Rcf. 



12|. The low-bias region corre- 



sponds to quasi-resonant direct transitions between the 
two Neel states. The intermediate-bias region is charac- 
terized by the threshold of chain excitations. Beyond this 
threshold, tunneling electrons induce an indirect process 
mediated by spin wave excitations. A second threshold 
defines the high-bias region where domain-wall excita- 
tions dominate the switching process. Correlation, i.e. 
mixing of spin configurations in the chain, is the driving 
force of the three Neel switching processes. However, cor- 
relation acts differently in the three, resulting in different 
dependences on the chain length of the three process and 
leading to the dominance of the process involving domain 
wall formation for long chains. 

The mechanisms described in this Letter are very dif- 
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ferent from the more usual local spin-flip mechanisms at 
play in, e.g., magnon excitation in ferromagnetic chains. 
Instead, the present mechanisms should be very general 
and operational in many systems with strong correla- 
tions, such as frustrated systems. 
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